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ABSTRACT 

We report on our serendipitous pre-discovery detection and detailed follow-up of the broad-lined 
Type Ic supernova SN 2010ay at z « 0.067 imaged by the Pan-STARRSl 37 t survey just ~ 4 days 
after explosion. Combining our photometric observations with those available in the literature, we 
estimate the explosion date and the peak luminosity of the SN, Mr « —20.2 mag, significantly 
brighter than known GRB-SNe and one of the most luminous SNe Ibc ever discovered. We measure 
the photospheric expansion velocity of the explosion from our spectroscopic follow-up observations, 
v p h ~19.2xl0 3 km s _1 at ~ 40 days after explosion. In comparison with other broad-lined SNe, the 
characteristic velocity of SN 2010ay is 2 — 5x higher and similar to the measurements for GRB-SNe 
at comparable epochs. Moreover the velocity declines two times slower than other SNe Ic-BL and 
GRB-SNe. Assuming that the optical emission is powered by radioactive decay, the peak magnitude 
implies the synthesis of an unusually large mass of 56 Ni, Mni = 0.9+q1 Mq. Our modeling of the 
light-curve points to a total ejecta mass, M e j ~ 4.7 Mq, and total kinetic energy, Er, 5 i ~ 11. Thus 
the ratio of Mn i to M e j is at least twice as large for SN 2010ay than in GRB-SNe and may indicate an 
additional energy reservoir. We also measure the metallicity (log(0/H) + 12 = 8.19) of the explosion 
site within the host galaxy using a high S/N optical spectrum. Our abundance measurement places 
this SN in the low-metallicity regime populated by GRB-SNe, and ~ 0.2(0. 5) dex lower than that 
typically measured for the host environments of normal (broad-lined) Ic supernovae. Despite striking 
similarities to the recent GRB-SN 100316D/2010bh, we show that gamma-ray observations rule out an 
associated GRB with E 1 < 6 x 10 48 erg (25-150 keV). Similarly, our deep radio follow-up observations 
with the Expanded Very Large Array rule out relativistic ejecta with energy, E > 10 48 erg. These 
observations challenge the importance of progenitor metallicity for the production of a GRB, and 
suggest that other parameters also play a key role. 

Subject headings: SurveystPan-STARRSl — gamma-rays: bursts — supernovae: individual (2010ay) 
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1. INTRODUCTION 

Recent observations have shown that long-duration 
gamma-ray bursts are accompanied by Type Ic super- 
novae (SNe) with broad absorption features (hereafter, 
“broad- lined” (BL)), i ndicative of high photosp heric ex- 
pansion velocities (see iWooslev fe Blo om 200(f for a re- 
view). This GRB-SN connection is popu larly expla ined 
by t he favored “collapsar model” dMacFadven et ahl 
l200l in which the gravitational collapse of massive 
(M > 20 Mg) progenitor star gives birth to a central 
engine - a rapidly rotating and accreting compact ob- 
ject - that powers a relativistic outflow. At the same 
time, not all SNe Ic-BL show evidence for a central en- 
gine. Radio observations constrain the fraction of SNe 
Ic-BL harboring relativistic outflow s to be less than a 
third (jSoderberg et al.lF2006al l2(7Toh . 

The physical parameter(s) that distinguish the pro- 
genitors of GRB-associated SNe from other SNe Ic-BL 
remains debated, while theoretical considerations indi- 
cate that progenitor m ctallicity may play a key role 
(IWooslev fc Hegerll2006D . In the collapsar scenario, mas- 
sive progenitor stars with metallicity above a thresh- 
old, Z > 0.3 Zq, lose angular momentum to metal 
line-driven winds, preventing the formation of a rapidly 
rotating compact remnant, and in turn, a relativistic 
outflow. At the same time, the hydrogen-free spec- 
tra of SNe Ic-BL indicate that their stellar envelopes 
have been stripped prior to explosion, requiring higher 
metallic ities (e.g., Z ~ Z& ) if due to radiation driven- 
winds (jWooslcv et alJfl995ll . Alternatively, short-period 
(~ 0.1 days) binary interaction may be invoked to spin up 
stars via tidal forces as well as cause mas s loss via Roche 
lobe overflow (jPodsiadlowski et~ahl 12004 iFrver fc Hegerl 
1200511 . However, even in the binary scenario, GRB for- 
mation is predicted to occur at higher rates in lower- 
metallicity environments, where the radius and mas s loss 
rates of stars should be smaller (llzzard et al.ll2004j) . 

Observationally, most GR B-SNe are discovered within 
dwarf star- forming galaxies (jFruchter et ali [2006 1 char- 
acterized by sub-sola r metallicities, Z < 0.5 Z Q 

(|Levesaue et all l2010al) . This has been interpreted as 
observational support for the metallicity-dependent col- 
lapsar model. Meanwhile, SNe Ic-BL without associated 
GRBs have historically be en found in more enriched en- 
vironments dModiaz et al.ll2008 1. motivating the sugges- 
tion of an observationally determined “ cut-off meta llic- 
ity” a bove which GRB-SNe do not form dKocevski et al.l 
I2009T) . However, this difference may be partly at- 
tributable to the different survey strategies: SNe have 
been found in large numbers by galaxy-targeted sur- 
veys biased towards more luminous (and therefore higher 
metallicity) environments, while GRB host galaxies are 
found in an untargeted manner through their gamma-ray 
emission. 

Against this backdrop of progress, recent observations 
have begun to call into question some aspects of this 
scenario. First, several long-duration GRBs have now 
been identified in solar or sup er-solar metallicity en - 
vironments (e.g., GRB 020819; ILevesaue et all l2010blh 
Similarly, the luminous radio emission seen from SN 
Ic-BL 2009bb pointed unequivocally to the production 
of co pious relativistic ejecta resembling a GRB after- 
glow (jSoderberg et al.l 120101 ) while the explosion envi- 


ronment was characterized by a sup er-solar metallicity, 
Z ~ 1 — 2 Zq (ILevesaue et al.ll2C)10a) . Together with the 
growing lack of evidence for massive progenitor stars for 
SNe Ic in pre -explosion Hubble Space Telescope images 
( Smarttl l2009( l . a lower mass (M ~ 10 — 20 Mq) binary 
progenitor system model (with a gentler metalli city de- 
pend ence) is gaining increasing popularity (lYoon et al.l 
I2010T) . Multi- wavelength studies of SNe Ic-BL in metal- 
poor environments may shed further light on the role of 
metallicity on nature of the progenitor and the explosion 
properties, including the production of relativistic ejecta. 

Fortunately, with the recent advent of wide-held opti- 
cal transient sur veys (e.g., C atalina Real-time Transient 
Survey ('CRTSh lDrake et al.ll2009L Panchromatic Survey 
Telescope and Ra pid Response Syste m; (Pan-STARRSl, 
abbreviat ed PS 1) Kais cret, al.l I2002L Palomar Transient 
Factory; lLaw et aLfl2nn9H SNe Ic-BL are now being 
discovered in metal-poor environmen ts. Z ~ 0.5 
(jArcavi et alJl20lfl iModiaz et al .1120 lTtl thanks to an un- 
biased search technique. In this paper, we present pre- 
discovery Pan-STARRSl imaging and multi- wavelength 
follow-up obser vations for the SN Ic-BL 2010ay discov- 
ered by CRTS (Drake et alJl2010 1. In (J21 we report our 
optical (Pan-STARRSl, Gemini, William Herschel Tele- 
scope) and radio (Expanded Very Large Array; EVLA) 
observations. In we model the optical light-curve and 
analyze the spectra to derive the explosion properties of 
SN 2010ay. In jj5j we use our observations of SN 2010ay 
with the EVLA to place strict limits on the presence 
of relativistic outflow. In (J5J we draw from gamma-ray 
satellite coverage to rule out a detected gamma-ray burst 
in association with SN 2010ay. In (23 we derive the ex- 
plosion site metallicity and find it to be significantly sub- 
solar and typical of most GRB-SNe host environments. 
In fJTJ we discuss the implications of these findings in the 
context of the explosion and progenitor properties and 
conclude in jJ5] 


2. OBSERVATIONS 
2.1. Discovery by CRTS 

SN 2010ay was discov ered b y the Cata lina Real- 
time Transient Su rvey (IDrake et al.l 120091 ) on 2010 
March 17.38 UT (IDrake et al.1 1201M and designated 
CSS100317:123527+270403, with an unfiltered magni- 
tude of m « 17.5 mag and located < 1” of the cen- 
ter of a compact galaxy, SDSS J123527. 19+270402.7 at 
z = 0.067 (Table [T]). We adopt a distance, Dl = 297.9 
Mpc, to the host galaxj£3, and note that the Galac- 
tic extinction toward this galaxy is E(B — V) = 0.017 
(jSchlegel et al.l fT998h . Pre-discovery unfiltered images 
from CRTS revealed an earlier detection of the SN on 
Mar 5.45 UT at m ss 18.2 mag a nd a non-detection from 
Feb 17.45 UT atm> 18.3 mag (jDrake et al.ll2010l) . 

A spectrum obtained on Mar 22 UT revealed the SN to 
be of Type Ic with broad features, similar to the GRB- 
associ ated SN 1998bw spectru m obtained near maximum 
light dFilippcnko et all 120101) . This classification was 
confirmed bv IPrietol ( 20101) who additionally reported 
photometry for the SN (see Table H). After numerically 
subtracting the host galaxy emission, they estimate an 
unusually bright absolute magnitude of V ~ —19.4 mag. 


19 We assume Hq = 71 km s 1 Mpc 1 , Pa = 0.73, Pm 


= 0.27 
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TABLE 1 

SN 2010AY HOST GALAXY SDSS 
J123527. 19+270402. 7 


Parameter 

Value 

RA 

12 h 35 m 27 s .2 (J2000) 

DEC 

+27° 04'03" (J2000) 

Redshift (z) 

0.0671 

Petrosian radius 

1.355" 

Photometry a 

u' 

19.56 ± 0.03 mag 

g' 

19.02 ± 0.01 mag 

r’ 

19.02 ± 0.01 mag 

i' 

18.69 ± 0.01 mag 

z ' 

18.87 ± 0.04 mag 

u 

19.50 ± 0.06 mag 

B 

19.02 ± 0.05 mag 

V 

19.13 + 0.05 mag 

R 

18.94 ± 0.05 mag 

I 

18.90 ± 0.06 mag 

Extinction 

E(B-V) MW b 

0.017 mag 

E(B-V) hoat c 

0.2 mag 

Note. — SDSS 

host galaxy properties 


and ugriz photometry. 

a Model magnitudes from SDSS DR8 
JAihara et al.| [20111) . Host galaxy photom- 
etry has not been corrected for extinction. 

U BV RI photometry has been converted 
from the SDSS ugriz measurements us- 
ing the transformation of lBlanton Sz Rowe id 
( 2002 ). 

k The Milky Way extinction as determined 
by iSchlegel et al.l (19981) , assuming Ry = 

3.1. 

c The host galaxy extinction determined 
from the SDSS spectrum centered on the 
galaxy nucleus, via the Balmer decrement 
as described in <j2.3. ll 

2.2. Pre-Discovery Detection With Pan-STARRSl 3n 

The field of SN 2010ay was serendipitously observed 
with the PS1 37 t survey in the weeks preceding its 
discovery. Pan-STARRSl is a wide-held imaging sys- 
tem at Ha leakala, Hawaii dedicated to survey obser- 
vation^ dKaisc r et al.l I2002D . The PS1 optical design 
(iHodann et al.l I2004T) uses a 1.8 meter diameter //4.4 
primary mirror, and a 0.9 m secondary. The telescope il- 
luminat es a diameter of 3.3 deg rees. The Pan-STARRSl 
imager ( Tonrv fc Onakal 120091 1 comprises a total of 60 
4800 x 4800 pixel detectors, with 10 ^m pixels that sub- 
tend 0.258 arcsec. The PS1 observations are obtained 
through a set of five broadband filters designated as gp±, 
rpi, ipi, zpi, and ypi. These hlters are simil ar to those 
used in previous surveys, such as SDSS (jFukugita et al.1 
119961 1. However, The gpi hlter extends 20 nm redward 
of g', the zpi hlter is cut o ff at 930 nm and S PSS has 
no corresponding ypi hlter dStubbs et alJl2010h . 

The held of SN 2010ay was observed on 2010 Febru- 
ary 21st (rpi-band) and February 25th (jpi-band, Fig- 
ure [TJ. On each night, four exposures were collected fol- 
lowing the strategy of the PSl-3pi survey (Chambers et 
al., in preparation). Following the CRTS discovery and 
announcement of SN 2010ay, we geometrically registered 
SDSS pre-explosion images to the PS1 images and per- 
formed digital image subtraction using the ISIS package 


TABLE 2 

SN 2010AY LIGHT-CURVE PHOTOMETRY 


MJD 

tpeak 

Filter 

m b 

Mr* 

Source 

55244.4 

-29 


< 18.3 

< -19.2 + 0.2 

d 

55248.6 

-25 

rpi 

< 22.0 + 0.1 

< -16.0 + 0.1 

e 

55252.2 

-21 

Zpi 

21.1 + 0.3 

-16.8 + 0.3 

e 

55260.4 

-13 


18.2 

-19.3 + 0.2 

d 

55272.4 

-1 


17.5 

-20.4 + 0.2 

d 

55277.2 

4 

B 

18.39 + 0.05 


f 

55277.2 

4 

V 

17.61 + 0.05 


f 

55277.2 

4 

R 

17.44 + 0.05 

-20.22 ± 0.07 

f 

55287.0 

14 

R 

18.2 + 0.2 

-19.2 + 0.2 

g 

55297.0 

24 

g 

18.9 + 0.1 


h 

55297.0 

24 

r 

18.3 + 0.1 

-19.0 + 0.1 

h 

55297.0 

24 

i 

18.0 + 0.1 


h 

55645.0 

372 

Zpi 

< 22.2 ± 0.2 

< -16.22 + 0.1 

e 

55649.0 

376 

r P i 

< 21.9 + 0.1 

< -15.7 + 0.1 

e 


a Time since peak in days, relative to the fitted value: 2010 March 18.1±1.9. 
k The measured apparent magnitude of the source, in the filter noted and 
without extinction correction. For the Pan-STARRSl photometry, a tem- 
plate image was subtracted; for the other points, the host galaxy flux has 
not been subtracted numerically. 

c The absolute R magnitude of the SN. Filter conversion, host flux sub- 
traction, and extinction correction have been performed as described in 
§exp:lcmodel. 

iDrake et al.l (20101) , unfiltered (synthetic V-band). 
e Photometry from the Pan-STARRSl 37 r survey. 
f IPrietol 112010) . 

® Synthetic photometry obtained from our WHT spectrum as described in 

pm 

Photometry from our Gemini/GMOS observations described in £12.21 


(! Alardl l2000f) . No residual flux was found in the differ- 
ence rpi-band image from 2010 February 21 with an up- 
per limit of r > 22.0 mag. However, we detect residual 
flux at the position of SN 2010ay in the *pi-band resid- 
ual image from 2010 February 25 with a magnitude of 
i = 21.1 ± 0.3 mag. 

The field was again observed in the ipiband on 2011 
March 25 and the rpi-band on 2011 March 29, but no 
residual flux was detected in the subtractions at the SN 
position to limits of i' > 22.2 and r' > 21.9. 

In Tableland Figure[2l we compile photometry from 
the PS1 detections, our optical observations, and the cir- 
culars to construct a light-curve for SN 2010ay. 


2.3. Optical observations 

We obtained an optical spectrum (~ 3000—11000 A) of 
SN 2010ay on April 1 UT, from the ISIS blue arm instru- 
ment of the 4.2 m William Herschel Telescope (WHT) 
at the Roque de Los Muchachos Observatory. The spec- 
trum was taken at the parallactic angle and the expo- 
sure time was 1800 sec. We obtained a second 1800 sec 
optical spectrum (~ 3600 — 9600 A) using the Gemini 
Multi-Object Spectrograph (GMOS) on the 8.1 m Gem- 
ini North telescope on 2010 April 11.4 UT. We employed 
standard two-dimensional long-slit image reduction and 
spectral extraction routines in IRAhPl We do not apply 
a correction for atmospheric differential refraction, be- 
cause the displacement should be < 0.5' at the airmass 
of the observations, « 1.0. 

20 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. - Images illustrat ing; the Pan-STARRSl pre-discovery detection of SN 2010ay. a) Pre-explosion i-band image from the SDSS 
DR7 llAbazaiian et aLll2009n . observed 2004 December 21, b) ipj-band image from the 3 tt survey of PS1, observed 2010 February 25, c) 
the difference of the SDSS i' and PS1 ipi images. The transient emission can be seen in the NE corner of the frame. 



and blueshift of the feature near 6355 A are similar for 
SN 2010ay and SN 2010bh, and are broader and more 
blueshifted than in SN 1998bw at a comparable epoch. 
We discuss the comparison between these two SN fur- 
ther in H3.3l and > 16.21 

Additionally, we obtained 60 s gri - band images of 
SN 2010ay on 2010 April 11.4 UT' using GMOS. The 
data were reduced using the gemini package in IRAF, 
and photometry was performed using the standard 
GMOS zero-pointO We measure that [g,r, i] = 
[18.90, 18.32, 18.04] ±0.1 mag. 

Imaging photometry is not available at the epoch of our 
WHT spectrum. We perform synthetic photometry on 
the spectrum to extract the flux at the central frequency 
of the R-band (6527 A) and find R = 18.2 ±0.2 mag. We 
then subtract the host galaxy flux numerically. 


Fig. 2. — The optical R-band light curve of SN 2010ay, as com- 
piled in Table[2]from CBET 2224 (red squares), the PS1 3 it s urvey 
(blue), our Gemini images (black, and synthetic photometry 1 32.31) 
based on our WHT spectrum (gold). Triangles denote upper limits. 
The thick dashed line represents the lumi nosit y of an expanding 
fireball fit to our early-time photometry ( 33.11). The thin da shed 
line is the SN Ibc light-curve template of ID rout et all (120101 1 and 
the gray field represents the standard deviation among its con- 
stituent photometry. The template is stretched by (1 + z) = 1.067 
with the best fit parameters fR pea k = Marchl8±2 (MJD 55273±2), 
AAtpeak = —20.2 zb 0.2 mag, and reduced x 2 = 0-5. The solid line is 
the light-curve of the SN-GRB 1998bw ( Galama et al.lfl998l '). red- 
shifted to match SN 2010ay. The vertical dotted lines mark the 
epochs of our Gemini and WHT spectroscopy. 


In both our Gemini and WHT spectra, broad absorp- 
tion features associated with the SN are clearly detected 
in addition to narrow emission lines typical of star form- 
ing galaxies. We distinguish the host galaxy emission 
from the continuum dominated by the highly broadened 
SN emission by subtracting a high-order spline fit to 
the continuum. Both SN and host galaxy spectral com- 
ponents are shown in Figure [3] As illustrated in the 
Figure, the broad, highly-blended spectral features of 
SN 2010ay resemble those of the Type Ic-BL SN 2010bh 
(associated with GRB1 00316D) at a similar similar epoch 
(jChornock et al.l I2010D . In particular, the broadening 


2.3.1. Host galaxy features 


We measure fluxes of the narrow emission lines from 
the host galaxy in our Gemini spectrum, as reported in 
Table (3| We fit a Gaussian profile to each narrow emis- 
sion line; for nearby lines such as [N II] and Ha, we fit 
multiple Gaussians simultaneously. We model the local 


continuum with a linear fit to 20 A regions on either side 
of each line. 

The host galaxy is significantly reddened as evidenced 
by the flux ratio of Hq, to H^, « 3.44. We infer E(B- 
V)= 0.2 mag {Ay = 0.6 mag), as measured from the 
Balmer decrement in our Gemini spectrum, assuming 
Rv = 3.1, Case B recombination ( Osterbrock fc Ferlandl 
I2006D . and the reddening law of Cardelli et al.l ( 19891) . 
This is similar to the value derived from the SDSS DR8 
nuclear fiber spectroscopy line fluxes for t he host galaxy 
(E[B-V)= 0.2 mag). The value reported in IModiaz et al.l 
(|2010f ) was also similar: E(B-V)= 0.3 mag. Furthermore, 
we note th at the color ( B — V = 0.78 ± 0.07 mag) as re- 
ported by iPrietol (120101) at 4 d after R— band peak is 
significantly redde r than SN Ibc color curve templates 
iDrout et all (120101) . further supporting a non- negligible 
host galaxy extinction. 


21 http:/ /www. gemini. edu/?q=node/10445 
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Fig. 3. — Optical spectra of SN 2010ay from Gemini (24 days after R-band peak) and the WHT (14 days). The spectrum is plotted 
decomposed into SN (above, with narrow lines clipped) and host galaxy (below, from Gemini, with spline- fit subtracted) components for 
clarity. The spectrum of SN 2010bh fromjChomock et_alJ (l2010l l is given in black for comparison, at 21.2 da ys after the GRB 100316D 
trigger (~ 10.0 days after R-band peak. lGanT‘et**ahn20Tlal l . The spectrum of SN 1998bw at +19 days from iPatat et ahl (120011 ') is also 
plotted. Both are transformed to the redshift of SN 2010ay. The SNe are shifted in flux for clarity. In the lower plot, relevant host galaxy 
emission lines are marked with a red line and labeled. 


2.4. Radio Observations 

We observed SN 2010ay with the EVLA dPerlev et alJ 
120091) on three epochs, 2010 March 26, 2010 April 29, 
and 2011 May 7. All EVLA observations were obtained 
with a bandwidth of 256 MHz centered at 4.9 GHz. We 
used calibrator J1221+2813 to monitor the phase and 
3C286 for flux calibration. Data were reduced using the 
standard packages of the Astronomical Image Processing 
System (AIPS). We do not detect a radio counterpart to 
SN 2010ay in these observations and place upper limits 
of F v <46, 42, 30 //Jy (3er) for each epoch respectively 
corresponding to upper limits on the spectral luminosity 
spanning L „ < (3.6 — 5.5) x 10 27 erg cm -2 s _1 . 

As shown in Figure [4j these limits are compa- 
rable to the peak luminos i ties o b served for ordi- 
nary SNe Ibc (j Berger et al.1 l2003al iSoderbergl 120071 : 


I Chevalier fc Franssonl 120061 iSoderberg et all 120101 and 
references within) and a factor of 10 * 2 to 10 3 less luminous 
than the radio afterglows associated with GR Bs 020903 , 
030329, and 031203 at ea rly epochs ( Berger et al.ll2003bt 
ISoderberg et al.l f2004allbl : iFrail et all 120051) . In compar- 
ison with the radio l uminosities observed for the rela- 
tivistic SNe 1998bw dKulkarni et ahl 119981 ) and 2009bb 

(|Soderberg et al.l 12010T) . SN2010av is a factor of > 10 
less luminous. The only relativistic explosion with de- 
tected radio emission below o ur EVLA limits is the w eak 
and fast fading XRF 060218 (ISoderberg et alJfeOOOal) . 

3. INITIAL CONSTRAINTS 
3.1. Light Curve Modeling 

We construct an R-band lightcurve for SN 2010ay us- 
ing the observations described in (12.31 (Table l2l). We con- 
vert the ipi band data points to the R — band assuming 
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Fig. 4. — EVLA upper limits for SN 2010ay (black arrows) are compared with the observed radio light-curves for ordinary SNe Ibc (grey; 
Soderberg 2003 and references withi n) and the radio aftergl ows of all GRB-SNe within z < 0.25. SN 2010ay is a factor of 10 2 to 10 3 less 
luminous than XRF 020 903 (orange; ISoderberg et aLll2004aT) , GRB 03 0329 (blue; IBerger et~ al . 2003b; Frail et aLll2005|), and GRB 031203 
(ISoderberg et al.ll2004bl) . Relativistic, engine-driven SNe 1998bw (red; IKulkarni et al.lll998l ) and 2009bb (dark blue; ISoderberg et alj|20i0h 
are a factor of 10 more luminous than the SN 2010ay limits on a comparable timescale, while XRF 060218 lies a factor of a few below the 
limits. We constrain the radio counterpart to be no more luminous than XRF 060218 and comparable to the peak luminosities of ordinary 
SNe Ibc. 


TABLE 3 

Emission line fluxes measured for the 
HOST GALAXY OF SN 2010AY 


Emission Line 

Flux 


(10 -16 erg s -1 cm -2 ) 


O II] A3726, 3729 

52 ± 2 

H<5 

4.5 ±0.2 

Hy 

9.9 ±0.1 

O III] A4363 

0.7 ±0.1 

H/3 

24.5 ±0.1 

O III] A4959 

30.4 ±0.1 

O III] A5007 

90.1 ±0.2 

N IIJA6548 

2.2 ±0.1 

Ha 

84.5 ±0.1 

N II]A6584 

6.33 ± 0.06 

S II]A6717 

7.46 ± 0.06 

'S II]A6731 

5.59 ± 0.06 


Note. — All fluxes have been measured 
from our Gemini spectrum. No reddening cor- 
rection has been applied. There is an addi- 
tional systematic uncertainty in the flux mea- 
surements of 10% due to flux calibration. 


the unextincted i' — R color observed fo r the Type Ic- 
BL SN 1998bw at the appropriate epoch dGalama et al.l 


Il998f) . Phot ometry for the unfi ltered CRTS images was 
reported by [Drake et all (l2010fl after transformation to 
the synthetic V-band (A.J. Drake, private communica- 
tion); we therefore assume the V — R color of SN 1998bw 
at the appropriate epoch. For the Pan-STARRSl pho- 
tometry, the host galaxy flux was subtracted using a 
template image. For all other photometry, we have sub- 
tracted the flux of the host galaxy numerically assuming 
the magnitude reported in Table |T] A total (Galactic + 
host) reddening of E(B-V)= 0.2 mag has been assumed 
(see H2.3.1D . 

To estimate the explosion date of SN 2010ay we have 
fit an expanding fire ball model to the op tical light curve 
(Figure [5]), following iConlev et al.1 (12006T) . In this model, 
the luminosity increases as 

We derive an explosion date to of 2010 February 
21.2±1.5. Here we have assumed an index n = 2. This 
suggests that the PS1 37r ipiband pre-discovery detec- 
tion image of SN 2010ay was taken ~ 4 days after the 
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explosion. This early observation provides a valuable 
datapoint for estimating the explosion date and also for 
constraining the rise-time of the SN. 

In Figure [21 we compare the light curve of SN 20 10ay 
to the SN Ibc light-curve template of lDrout ct al. 1 201 01 
after stretching by a factor of (1 + z). The template pro- 
vides a reasonable fit to optical evolution of SN 2010ay. 
Fitting the template to our photometry we derive (re- 
duced x 2 = 0-5) a date of R-band peak of 2010 March 
18±2 UT (MJD 55273±2) and a R-band peak magnitude 
of Mr ss —19.7 mag before extinction correction. As dis- 
cussed in 112.3.11 based on the Balmer decrement observed 
for the host galaxy emission lines, we assume an extinc- 
tion of E(B — V) = 0.2 mag. Applying this correction, 
the peak absolute magnitude is Mr S3 —20.2 ± 0.2 mag. 
We note that this fitted value is ss 0.2 mag fainter than 
that estimated from the data point near peak. Here, the 
uncertainty is dominated by the template fitting. 

Regardless of the extinction correction, SN 2010ay 
is brighter tha n all the 25 SNe Ibc in the sample 
of IDrout et al.1 (120101) . except for SN 2007D (Mr « 
—20.65 mag, which was also significantly extincted: 
Ay ~ 1.0 mag). Assuming an intrinsic V-R col or 
of zero at peak (e.g. 1998bw: [Galama et al.l 119981 : 

iPatat et al.l l200lh . SN 2010ay is also brighter than any 
of the 22 GRB and XRF-p roducing SNe in the compila- 
tion of lGano et al.l ( 2011b ). all corrected for extinction, 
and is 1.5 standard deviations from the mean bright- 
ness. The peak magnitude is only ~ 1 mag below that 
of the Type Ic SN2007bi (Mr = -21.3 ± 0.1 mag), 
which [Gal- Yam ct al.1 ([2009 ) report as a candidate pair- 
instability supernova. 

3.2. Large Nickel Mass for SN 2010ay 

We use the available photometry for SN 2010ay dis- 
cussed above to derive the mass of 56 Ni required to power 
the optical light-curve under the assumption that the 
emission is powered by radioactive decay. Using the 
relat ion between M^i and Mr found by IDrout et al.l 
(12010D . log(MNi) = (—0.41 Mr — 8.3) M 0 , we esti- 
mate that SN 2010ay synthesized a nickel mass of Mn i = 
0.9+ob Mq. We have estimated the uncertainty in the 
Mm estimate by propagation of the uncertainty in the 
template fitted peak magnitude — systematic uncer- 
tanties are not included. If we instead adopt the most 
luminous individual data point in the light curve as the 
peak value, rather than the smaller peak value from tem- 
plate fitting, we find Mn; rs 1.2 Mq. 

The Mm estimate for SN 2010ay is larger than that 
of all but one (S N 2007D) of the 25 SNe Ibc in the 
IDrout et~ah ( 1201 01 compilation and significantly larger 
than the estimate for G RB-SN 2010bh, Mm = 0.12 ± 
O.OIMq, as estimated bv I Cano et al.l (1201 lal) . On the 
other hand, Mm of SN 2010ay is at le ast 3 x smaller tha n 
for SN 2007bi (M m « 3.5 - 4.5 M^. lYoung eIaP[20l?l . 
A pair instability supe rnova should produce Mm > 3 M 0 
(ICal-Yam et a,lJl2009l) . 

3.3. Unusually high photospheric expansion velocity 

As illustrated in Figure [3] the broad, highly-blended 
spectral features of SN 2010ay at the time of the WHT 
observations (14 days after R-band peak, see < 13.11) re- 
semble those of SN 2010bh at a similar time (10.0 d after 


peak. IChornock et al.ll2010h . In particular, the blueshift 
of the feature near 6355 A is larger than in SN 1998bw 
and more similar to SN 2010bh. This feat ure is com- 
monly associated with Si II A6355A (e.g. IPatat et~ahl 
1200 11 1. However, this feature has two clearly-detectable 
absorption minima in SN 2010ay, but not in SN 2010bh. 
This could be due to increased blending in SN 2010bh 
or the absence of contaminating lines. The red ends of 
the SN 2010ay and SN 2010bh spectra (rest wavelength 
> 7500 A) have similar P-Cygni features, but the emis- 
sion and absorption components in SN 2010bh are each 
blueshifted by ~ 200 A m o re th an in the spectrum of 
SN 2010ay. IChornock et al.l ([2010D attribute this feature 
to the Ca II NIR triplet, with a gf - weighted line cen- 
troid of 8479 A, and find a velocity that is high, but 
consistent with the early-time velocity of Si II A6355A 
(30- 35 x 10 3 km s" 1 ). 

We estimate the photospheric expansion velocity (v p h) 
from the minimum of the Si II A6355A absorption 
feature. We smooth the spectru m using an inverse- 
variance-weighted Gaussian filter dBlondin et all 120061 . 
with dX/X = 0.005) and measure the minimum position 
of the redmost component of the absorption profile. The 
blue component of the absorption profile shifts blueward 
over time, suggesting that it is produced by a combina- 
tion of ions separated by several 10 3 km s _1 , such as 
He I A5876A and Na I D, whose relative optical depth 
changes with time. 

The photospheric velocity inferred from the 
Si II A6355A feature is ~ 2x faster than that of SN 
1998bw at similar times, and more simi lar to that of SN 
2010bh (Figure 3). For SN 1998bw. IPatat et all ([20011 1 
measured ~ 10 x 10 3 km s -1 at +13 days. For SN 2010bh, 
IChornock et al.l (|2010l) measured v p h ~ 26 x I Q 3 km s _1 
at +10.0 days after explosion. iPrietol (IMol) re- 
ported a velocity of v p h ~22.6xl0 3 km s _1 from the 
Si II A6355A feature in a spectrum of SN 2010ay taken 
at +4 days. From our [WHT,Gemeni] spectra taken 
[+14, +24] days after R-band peak (see 1 13.11) . we 
estimate v p h ~ [19.2,18.3] x 10 3 km s _1 . 

In addition to the broadening of the spectral fea- 
tures and the blueshift of the Si II A6355A line, ad- 
ditional lines of evidence suggest a high photospheric 
expansion velocity for SN 2010ay. We measure v p h ~ 
[21.7, 20.1] x 10 3 km s _1 from the Ca II NIR triplet on the 
smoothed [WHT, Gemini] spectra, relative to a line cen- 
ter at 8479 A. This is within a few 10 3 km s _1 of the v p h 
we measure from Si II A6355A at these epochs. Further- 
more, we do not detect the broad emission “bump” near 
4500 A in either of our spectra. This feature was also 
absent in SN 2010bh, but was identified in SN 2003dh, 
SN2006aj, and seve ral Ic-BLs not as s ociate d with GRBs 
and normal SNe Ic; IChornock et al.l ( 2010 9 suggest that 
the absence of this feature indicates a high expansion ve- 
locity if it is due to blending of the iron lines to the blue 
and red of 4500 A. 

We compare the late-time expansion velocity of 
SN 2010ay to other SNe Ic-BL and GRB-SNe with de- 
tailed time-dependent velocity measurements from the 
literature in Figure[5] In this figure, we also fit power-law 
gradients to the time-evolution of the velocity of these 
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TABLE 4 

Velocity Evolution of 
SNe Ic-BL 


SN vf h 


SNe Ic-BL 


1997ef 

6 

-0.8 

2002ap 

4 

-1.1 

2003jd 

10 

-0.5 

2007bg 

7 

-0.2 

2007ru 

10 

-0.5 

2010ay 

22 

-0.4 

Engine-driven SNe Ic-BL 

1998bw 

10 

-0.9 

2003dh 

12 

-0.9 

20031w 

10 

-0.8 

2006aj 

15 

-0.3 

2009bb 

11 

-0.8 

2010bh 

24 

-0.2 


Note. — To the velocity 
measurements for each SN, 
we have fit a power law of 
the form v ph = v p ° h (t/ 30)“, 
where t is the time since ex- 
plosion in days and v 3 ^ is the 
velocity at 30 days in units of 
10 3 km s -1 . The parameter 
o; then represents the expo- 
nential velocity gradient and 
u 3 ^ is a characteristic veloc- 
ity. These power law fits are 
illustrated in Figure [5] 


SNe with the form v p h = u 3 °(t/ 30)“, where t is the time 
since explosion in days and v 3 ° is the velocity at 30 days 
in units of 10 3 km s _1 . These parameters are listed in 
Table |31 Figure [5] illustrates that most SNe are well de- 
scribed by a single power lawJ3 Contamination from 
different ions or detached features will add uncertainty 
to velocities measured from the Si II A6355A feature. 

SN 2010ay and 2010bh share high characteristic ve- 
locities at 30 days after explosion and velocity gradi- 
ents that are slow relative to other broad-lined Ic SNe 
with and without associated GRBs. For SN 2010ay, 
u 3 ° = 22 is 2 — 5x larger than for other SNe Ic-BL 
without associated GRBs (6 for 1997ef, 4 for 2002ap, 
10 for 2003jd, 7 for 2007bg, and 10 for 2007ru) and is 
similar to the GRB-SN 2010bh (u 3 ° = 24). No other 
GRB-SN or SN Ic-BL has vf h > 15. The SNe Ic-BL 
and GRB-SNe with the most shallow velocity gradients 
among these twelve objects have a < —0.5; they are 
SN 2006aj ( a = -0.3), SN 2007bg (a = -0.2), SN 2010ay 
(i a = — 0.4), and SN 2010bh ( a = —0.2). Two of these 
objects are GRB-SNe, one is SN 2007bg (whose small 
velocity at early times distinguishes it from other SNe 
Ic-BL), and the other is SN 2010ay. The velocity of SNe 
2010ay and 2010bh, respectively, declines about 2 and 
4x more slowly than the other SNe Ic-BL (mean and 
standard deviation: a = —0.8 ± 0.3) and about 1.5 and 
3x more slowly the GRB-SNe ( a = —0.6 ± 0.3). 

Given the high peak luminosity of SN 2010ay ( H3.1D . we 

22 A break at ~ 2 X 10 4 km s _1 appears to exist for SN 1998bw 
at ~ 16 days. 


also consider the velocity of the candidate pair-instability 
SN 2007bi. Velocity measurements for SN 2007bi are only 
available at late times (> 50 d. lYoung et al.lf2010h . Fit- 
ting to these late-time Si II A6355A velocity measure- 
ments, we find that SN 2007bi has a characteristic veloc- 
ity ~ 3x smaller than 2010ay: u 3 ° = 8 and the late-time 
velocity gradient is ~ 4x slower: a = —0.1. 


3.4. Ejecta Mass and Energy 

We use the scaling relations provided bv IDrout et al] 
(|2010t) . based on the original formali sm of lArne ttl (119821) 
and modified bv I Valenti et a l. (2008), to derive the total 
mass of the ejecta and the kinetic energy. 


M e j = 0.8 



2 


Vph 

10,000 km s- 1 


M 0 (2) 


E k = 0.5 



2 


Vph 
10, 000 km s 


-l 


3 

x 10 51 erg (3) 


We assume the fitt ed p eak magnitude for SN 2010ay 
(Mr = —20.2 ± 0.2, H3.ll) . the photospheric expansion 
velocity we measure from our WHT spectrum at 14 days 
after maximum light ( v p h = 19.2 x 10 3 km s -1 , if 1.311 . 
and a characteristic time (light-curve width) consistent 
with the data and the mean value from the iDrout et al.l 
(120101) sample of SNe Ic-BL (r c = 14 d). 

Using these values, the total mass ejected was M e j ps 
4.7 M 0 , and the total kinetic energy of the explosion 
was Er ps 10.8 x 10 51 ergs. Hereafter we refer to the 
definition Ek , 51 = E k / 10 51 ergs. 

The systematic uncertainties associated with this mod- 
eling dominate the statistical uncertainties. In particu- 
lar, the models rely on the assumptions of homologous 
expansion, spherical symmetry, all 56 Ni centralized at the 
center of the ejecta, optically thick ejecta and constant 
opacity. 

We note that an earlier measurement of the photo- 
spheric velocity is preferred for optical modeling. Since 
we have argued that SN 2010ay and SN 2010bh have sim- 
ilar characteristic velocities (u 3 °), if we instead adopt 
a higher velo city of 25 , 000 km s _1 as measured for 
SN 2010bh by I Cano et all (l2011al) . we estimate M e j pa 
6.1 M 0 and Ek , 51 ~ 23.9 for SN 2010ay. 

The M e j and Ek ,5 1 of SN 201 0ay are consisten t with 
the mean for SNe Ic-BL in the IDrout et alJ ( 2010 ) sam- 
ple (4.7+Lg M 0 and 1 1 A 4 , respectively), because the au- 
thors assumed a velocity ( v p h = 2 x 10 4 km s _1 ) similar 
to the late-time velocity we measure. For comparison, 
SN 2010bh had a total ejecta mas s of ~ 2 M w an d a 
total kinetic energy of Ek , 51 ~ 13 dCano et al.l 1201 la ). 

The ratio of Ni to total ejecta mass is ~ 0.2 for SN 
2010ay, significantly higher than the values typical of SNe 
Ic-BL and GRB-SNe. For comparison, the ratio is just 
~ 0.05 for SN 2010bh. Adopting th e values deri v ed from 
bolometric light curve modeling bv iCano et al.l ( 2011a ). 
the Mni and MNi/M e j ratios for other GRB-SNe are: ~ 
0.5 M 0 and ~ 0.06-0.22 (1998bw), ~ 0.4 M 0 and - 0.08 
(2003dh), ~ 0.15 M 0 and - 0.07 - 0.1 (2006aj), and ~ 
0.2 M 0 and ~ 0.06 (2009bb). This ratio for SN 201 0ay is 
larger than that of all but 4 of the 25 SNe of lDrout et al.l 
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Fig. 5. — A comparison of the time-evolving photospheric expansion velocity of SN 2010ay and other SNe Ic-BL (red) and GRB-SNe (blue) 
from the literature. Typical uncertainties in velocity estimates are 10%. For each SN, we fit a power law of the form v p h = vb)) (f/30) a , 
where t is the time in days, v is the velocity at 30 days since explosion (dashed vertical line), and a is the velocity gradient. The velocities 
for SNe 1997ef, 2003dh, 20031w are from lMazzali et al.l 12000 ). as determined by spectral modeling. The velocities for all other SNe are 
measured from t he Si II A6355A feature as follows: SN 2007ru are fromlSahu et all (120091) ; SNe 1998bw, 2006aj, and 2010bh are from 
IChornoc k et al.l 1201011. fr om spectra in references therein; SN 2007bg are from lYoung et alj 1201Cf l: SNe 2002ap, 2009bb, and 2003jd are 
from IPignatiTTit all 120111 ) and references therein. 


IM3): the Type Ic SNe 2004ge ( M m /M ei ~ 0.4) ans 
2005eo (M Ni /M ej - 0.2), the Type lb SN2005hg (~ 0.4), 
and the Type Ic-BL SN2007D (~ 0.6). 

The large value of Mni we estimate for SN 2010ay raises 
the question of whether a process other than Ni decay 
may be powering its light-curve. An independent test 
of the physical process powering the light-curve is the 
decay rate of the late-time light curve which should be 
« 0.01 mag day -1 for SN e powered by radioact ive de- 
cay of 56 Co. For SN 2007bi. lGal-Yam et alJ (|2009D derive 
Mj vi = 3.5 Mq from the measured peak magnitude and 
find that the late-time light curve is consistent with the 
decay rate of 56 Co. While the Pan-STARRSl 37 t sur- 
vey also observed the field in 2011 March, the SN was 
not detected in our subtracted images and the limits are 
not constraining in the context of 56 Co decay (see 112.21 
and Table El). Another possible process is a radiation- 
dominated shock that emerges due to interaction with an 
opaque c ircumstellar medium, as has recently been pro- 
posed by Cheval ier fc Irwiiil (120111) for the ultra-luminous 
SNe SN 2006gy and SN 2010gx dPastorello et alJ 120101 : 
IQuimbv et al.l 120111) . However, while this class of 
ultra-luminous obje cts shares some spectr oscopic sim- 
ilarities to SNe Ic dPastorello et al.1 1201 0). they show 
peak luminosities ~ 4 — 100 x higher than SN 2010ay 


(jChomiuk et al.ll2011l : IQuimbv et al.ll2C)Tll) . 

4. CONSTRAINTS ON RELATIVISTIC EJECTA 

We use our EVLA upper limits for SN 2010ay spanning 
At « 29 — 433 days to constrain the properties of the 
shockwave and those of the local circumstellar environ- 
ment. The radio emission from SNe Ibc and GRBs is pro- 
duced by the dynamical intera ction of the fast est ejecta 
with the surrounding material dChcvalierlll982D . The ki- 
netic energy of the ejecta is converted, in part, to internal 
energy of the shocked material which itself is partitioned 
between relativistic electrons (e e ) and amplified magnetic 
fields (es). Following the breakout of the shockwave from 
the stellar surface, electrons in the environment of the 
explosion are shock-accelerated to relativistic velocities 
with Lorentz factor, q e and distributed in a power-law 
distribution characterized by -/V( 7 e ) oc 7 ~ p . Here, p char- 
acterizes the electron energy index. The particles gyrate 
in amplified magnetic fields and give rise to non-thermal 
synchrotron emission that peaks in the radio and mm- 
bands in the days to weeks following explosion with ob- 
served spectral index, F v oc _ At lower frequen- 

cies the emission is suppressed due to synchr otron se lf- 
absor ption which defines a spectral peak, v p dChevalied 

EUl). 
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The dynamics of the shockwave determine the evo- 
lution of the synchrotron spectrum, and in turn, the 
properties of the observed radio light-curves. In the 
case of SNe Ibc, there are three primary scenarios for 
the dynamical regime of the ejecta depending on the 
shock velocity, v = /3c (associated Lorentz factor, T): 
(i) non-relativistic ( v ss 0. 2c) free-expansio n as in the 
case of ordinary SNe Ibc (IChevalierl I1998D . (ii) a de- 
coupled and relativistic (T ~ 10) shell of ejecta that 
evolves accordi ng to the Blandf ord-McKee solution for 
several months (ISari ct a l.||1998|) before tra nsitioning to 
the Sedov- Tavlor regime dFrail et al.ll2000D . This is the 
standard scenario for typical GRBs. And (iii) a sub- 
energetic GRB with trans-relativistic velocity (/3T < 3) 
that bridges the free-expansion a nd Bl a ndford McKee 
dynamical regimes (e .g. SN 1998bw; rKulkarni et aljH998l 
ILi fe Chevalier! IT999I) . 

We consider our EVLA upper limits in the context 
of these three models below. For shock velocities of 
v > 0.2 c, e R ~ 0.1 is reaso nable (|Soderberg et all [20051 
IChevalier fc Fransso n 2006). We further assume equipar- 
tition, e e = e B = 0.1. We adopt a free expansion model 
for both the non-relativistic ordinary SN Ibc case and the 
su b-energetic, trans-relat ivistic GRB scenario. As shown 
by ILi fc Chevalierl (119991 ) a free-expansion model is still 
reasonable in the trans-relativistic regime (cases i and iii, 
see above). 


4 . 1 . Freely- expanding shockwave 

In the free-expansion scenario, a shock discontinuity 
separates the forward and reverse shocks, located at the 
outer edge of the stellar envelope. The bulk ejecta is in 
free expansion while the thin layer of p ost shock materi al 
is slightly decelerated, R oc t 0 9 (IChevalier fc Franssonl 
12001 . At a given frequency, the bell-shaped light-curves 
of the SN sync hrotron emission may be described as 
(IChevalierl 1 19981) 


K 


1.582 x L„ iP 


/ AA a 

i _ e -(A t/ tp )-(“+ b y 

\W 

■ 


( 4 ) 


where L vp is the flux density at the spectral peak at 
epoch, t p . Assuming an electron i ndex of p ss 3, consis- 
tent w ith radio spectra of SNe Ibc (IChevalier fo Franssonl 
I2006D . the exponents are a ~ 2.3 and b rs 1.3. The time 
averaged shockwave velocity is v « R/At where R is the 
shockwave radius defined as 


R 


2.9 x 10 16 




10 28 erg s -1 Hz -1 


9/19 


( 5 ) 


Here we make the assumption that the radio emitting 
region is half of the total volume enclosed by a spheri- 
cal blastwave. Next, we estimate the internal energy, E, 
of the radio emitting material from the post-shock mag- 
netic energy density, E ~ B 2 R? /\2e B where w e maintai n 
the as sumption of equipartition. As shown bv IChevalierl 
(I1998D . the amplified magnetic field at peak luminosity 
is also directly determined from the observed radio prop- 
erties, 


B 


/ e \ 

-4/19 

0.43 ( — ) 


\e B J 


x ( Up ) 

G. 

\5 GHz ) 


J v,p 


10 28 erg s- 1 Hz 


-l 


- 2/19 


( 6 ) 


Finally, the mass loss rate of the progenitor star, 
M, may be derived from the number density of emit- 
ting electrons. Here we normalize the wind profile ac- 
cording to p oc Ar~ 2 and A* = A/5 x 10 11 g cm^. 
This normalization of A* implies that an A* of 1 corre- 
sponds to typical Wolf-Rayet progenitor wind properties 
of M = 10 -5 Mq yr -1 and a progenitor wind velocity of 
v w = 10 3 km s _I . 



where we assume a shock compression factor of ~ 4 and 
an nucleon-to-electron ratio of two. 

We built a two-dimensional grid of fiducial radio light- 
curves according to Eqn. [I] in which we vary the pa- 
rameters L Vp and v p over a reasonable range of pa- 
rameters space, bounded by t p « [1,3000] days and 
F U}P « 0.04 — 1000 mJy. We identify the fiducial light- 
curves associated with a radio luminosity higher than 
the EVLA upper limits for SN 2010ay at each epoch as 
these are excluded by our observations. We extract the 
physical parameters associated with these excluded light- 
curves ( R , B , E, A*) to define the parameter space ex- 
cluded by our radio observations. The parameter space 
for v p and F v p are bounded by the respective values for 
which the model exceeds relativistic velocities, /? T ~ few. 

As shown in Figure [6l our deep EVLA limits enable us 
to rule out a scenario in which there is copious energy, 
E > 10 48 erg, coupled to a relativistic outflow, in this 
two-dimensional E — v parameter space. The excluded 
region includes GRB-SNe 1998bw and 060218 as well as 
the relativistic SN 2009bb. It does not exclude the stan- 
dard scenario in which a small percent of the energy is 
coupled to fast moving material within the homologous 
outflow, as is typically observed for ordinary SNe Ibc 
(E Rs lO 47 and v ~ 0.2c: iSoderberg et al.112010 ). 

Next we consider the effects of circumstellar density 
since lower mass loss rates produce fainter radio coun- 
terparts. As shown in Figure [7] the EVLA limits for 
SN 2010ay exclude the region of parameter space popu- 
lated by SNe 1998bw and 2009bb with mass loss rates 
of A* ~ 0.1, however, the low density environment of 
GRB 060218 lies outside of our excluded region due to 
its lower CSM density, A* ~ 0.01 that gives rise to a 
lower luminosity radio counterpart. 

4 . 2 . Relativistic Ejecta 

In the case of relativistic deceleration the ejecta are 
confined to a thin jet and are physically separated from 
the homologous SN component. Deceleration of the jet 
occurs on a timescale of At m (E^ /A*) years in a 
wind-stratified medium (IWaxmanl 120041) . On this same 
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Blastwave Velocity (T/S) 

Fig. 6. — The region of energy- velocity space ruled out (red) by 
our EVLA observations for on- axis ejecta under the assumption of 
a free-expansion model. 

timescale, any ejecta components that were originally 
off-axis spread sideways into the observer’s line-of-sight. 
While the early EVLA limits constrain the properties of 
the on-axis ejecta according to the free-expansion model 
described above, the late time EVLA upper limits con- 
strain any radio emission from a GRB jet originally 
pointed away from our line-of-sight. 

For this scenario, we a dopt the semi-analytic model 
of lSoderberg et all l]2006cD for off-axis GRB jets and de- 
rive the two-dimensional parameter space (energy and 
CSM density) that is excluded based on our EVLA up- 
per limit at At w 1.2 years. We note that this model 
accommodates the full transition from relativistic to non- 
relativistic evolution. We built a collection of model 
light-curves spanning parameter range, A* ps [0.01 — 100] 
and E ss [10 49 — 10 52 ] erg, maintaining the assumption of 
p = 3 and e e = cb = 0.1. Here we adopted a jet opening 
angle of 0j = 5° and an off-axis viewing angle of 9 oa . = 90° 
(the most conservative scenario). As shown in Figure 0 
we are able to exclude the parameter space associated 
with typical GRBs, i.e. E s=s 10 51 erg (beaming cor- 
rected) and A* « 1. GRBs with lower energies and den- 
sities are better constrained using the trans-relativistic 
formalism above. In conclusion, our radio follow-up of 
SN 2010ay reveals no evidence for a relativistic outflow 
similar to those observed in conjunction with the nearest 
GRB-SNe, however a weak afterglow like that seen from 
XRF 060218 cannot be excluded. 

5. CONSTRAINTS ON AN ASSOCIATED GRB 

Given the estimate of the explosion date we have de- 
rived (» we have searched for gamma-ray emission 
that may have been detected by satellites. No GRBs 
consistent with SN 2010ay were reported in the circulars 
of the Gamma-ray Coordinates Network Circulars, but 
it is possible that bursts were detected below the instru- 
ment triggering thresholds. 

We next consulted the sub-threshold bursts from the 
Swift Burst Alert Tele scope (BAT; iGehrels et al.1 120041 : 
iBarthelmv et al.ll2005fl detected within the 6 days sur- 
rounding the explosion date estimate. We find that 
no gamma-ray emission was detected within 0.5 degrees 



Fig. 7. — The region of energy-mass loss space ruled out (red) by 
our EVLA observations for on- axis ejecta under the assumption of 
a free-expansion model. 



Fig. 8. — Our EVLA observations at At « 1.2 years after explo- 
sion constrain the properties of a possibly associat ed of f-axis GRB 
jet. Using our semi-analytic model as described in 3m we assume 
partition fractions of e e = e b — 0.1, Qj = 5 degrees, p = 2.5, 
and an off-axis viewing angle of 6 oa = 90 degrees. We are able to 
exclude the region of Ek — A* parameter space (red) associated 
with typical GRBs, i.e. Ek = 10 51 erg and A* = 1 (dashed black 
lines). 

of the position of SN 2010ay by the BAT during this 
period. Given the sensitivity of the BAT, this cor- 
responds to an upper limit on the gamma-ray flux of 
~ 10 -8 erg s _1 cm -2 (15-150 keV). However, the field 
of SN 2010ay was in the field of view of the instrument 
for only 106 ksec during these 6 days, or ~ 20% of the 
duration. 

For complete temporal coverage, we have searched 
the records of the interplanetary network (IPN), which 
is sensitive to bursts with fluences down to ~ 6 x 
1 0~ 7 er g c m~ 2 ( 25-150 keV) (50% efficiency limit, 
iHurlev et all |2010j ), and observes the entire sky with 
a temporal duty cycle close to 100%. An undetected, 
sub-threshold burst should have a fluence below this 
limit. Between 2010 February 21 and 25, inclusive, a 
total of 12 bursts were detected by the spacecraft of the 
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IPN (Mars Odyssey, Konus-Wind, RHESSI, INTEGRAL 
(SPI-ACS), Swift-BAT, Suzaku, AGILE, MESSENGER, 
and Fermi (GBM)). Ten of them are confirmed bursts; 
they were observed by more than one instrument on 
one or more spacecraft, and could be localized. Two 
of them are unconfirmed bursts; they were observed by 
one instrument on one spacecraft (Suzaku). The total 
area of the localizations of the 10 confirmed bursts is 
~ 0.58 x 47 t sr. This implies that about 0.58 bursts can 
be expected to have positions that are consistent with 
any given point on the sky simply by chance (i.e. within 
the 3cr error region), and indeed none of the bursts in this 
sample has a position consistent with the SN position. 

These non-detections imply upper-limits to the 
gamma-ray energy (U 7 ) of a burst that may have been 
associated with SN 2010ay. The IPN non-detection in- 
dicates E-y < 6 x 10 48 erg (25-150 keV), while the BAT 
non-detection indicates that the peak energy of the burst 
was < lx 10 47 erg s _1 (15-150 keV) if the burst occurred 
while in the field of view of the instrument. 

By comparison to other GRB-SNe, it is unclear 
whether the gamma-ray emission from a GRB associated 
with SN 2010ay would have been detected by BAT or the 
IPN. The isotropic prompt energy release of long GRBs is 
typically ; so ~ 10 52 erg, however the prompt emission 
of the sub-energetic class of GRB-SNe can be sev eral or- 
ders of magnitude fainter dSoderberg et al.ir2006bf) . GRB 
980425/SN 1998bw had a peak gamma-ray luminosity of 
~ 5 x 10 46 erg s _1 (24-1820 keV iGalama et all 119981) . 
which is a factor of 2 below our BAT limit, and E lt ; so ~ 
5 x 10 47 erg dPian et all 12000 ,1. more than a factor of 
five below our limit. Neither satellite should have de- 
tected such a burst. In contrast, GRB 031203/SN 20031w 
had a peak gamma-ray luminosity of ~ 1 x 10 49 erg s _1 
(20-200 keV) and a total isotropic equivalent energy 
of i so = (4 ± 1) x 10 49 erg (iSazonov et al.1 120041) . 
about two orders of magnitude above the sensitivity of 
the BAT and twice the threshold of the IPN, respec- 
tively. GRB 030329/2003 dh was even more luminous, 
with E ly iso ~ 7 x 10 49 erg dHiorth ct all 120031 ) . A burst 
like GRB 031203 or 030329 should certainly have been 
detected by IPN, or the BAT if it occurred while the 
field of SN 2010ay was in the field of view of the in- 
strument. XRF 060218/SN2006aj, an extremely long- 
duration (At « 2000 s) event, had a peak luminosity 
observed by BAT of ~ 10 -8 erg cm -2 s _1 (15-150 keV), 
corresponding to a peak emission of « 2.4 x 10 46 erg s _1 
given the redshift of the burst (z=0.033), and a total 
isotropic equivalent energy of i so = (6.2 ± 0.3) x 
10 49 erg flCammna et al.ll2006l) . If such a burst was asso- 
ciated with SN 2010ay, its peak emission may have been 
a factor of four below the BAT sensitivity limit, although 
its total isotropic energy emission is an order of magni- 
tude larger than our IPN limit for SN 2010ay. Finally, the 
event whose host galaxy and supernova properties seem 
most similar to SN 2010ay, GRB 100316D/SN 2010bh 
had E 1 } i so > (5.9 ± 0.5) x 10 49 erg iStarling et alJ l)2011h 
— a full order of magnitude above our IPN limit. 

Another possibility is that prompt emission associated 
with SN 2010ay may have been too soft to be detected 
by the BAT or IPN. For example, the spectrum of XRF 
060218 rose to a peak at 0.3-10 keV at ~ 985 s after 
triggering, then softened significantly thereafter. Even 


though the the total emission of this burst is well above 
our IPN limit, it may have escaped detection if it was 
similarly soft. 

6. SUB-SOLAR HOST ENVIRONMENT METALLICITY 

We estimate the oxygen abundance of the host environ- 
ment of SN 2010ay from the strong nebular emission line 
fluxes measured from our Gemini spectrum (Table[3|). At 
the distance of the host galaxy, the 1' Gemini slit width 
corresponds to a physical size of 1.3 kpc. 

From the 03N2 diagnostic of lPettini fc Pavel (I2004T ) 
(PP04), we derive a metallicity of log(0/H) = 8.19, or 
Z ~ 0.3 Zq, ad opting the solar meta llicity log(O/H) 0 ± 
12 = 8.69 from lAsplund et ai~l (j2005t ). Using the N2 di- 
agnostic of PP04, we find log(0/H) + 12 = 8.26. Us- 
ing the abundance diagnostic, R 23 = log([0 II] A3727 ± 
[O II I1AA4959,5007)/H& we find log(0/H)+12 Z94 = 
8.49 ( Zaritskv et all I1994D a nd log(0/H)±12 KD 02 = 
8.51 (IKewlev fe Dopital 120021) . However, these R 23 - 
based estimates are more sensitive to flux-calibration 
and reddening-correction. Moreover, there is a well- 
known bi-valued relationship between R23 and oxy- 
gen abundance. The R23 value of the explosion site 
(0.872±0.003) places it near the turnover point, but 
we assume that it lies on the up per branch based 
on its [N II] /[O II] ratio, following IKewlev fc Ellisonl 
(120081 ). The metallicity values we derive using the 
PP04 and KD02 diagnostics are approximately equiv- 
alent giv en the offset that exists between these two di- 
agnostics lKewlev fe Ellisonl (120081) . These measurements 
are similar to the values reported bv lModiaz et all (|2010f ) 
(log(0/H)+12 [PP04,KD02]=[8.2,8.4]) for the SN2010ay 
host galaxy. The statistical errors in our strong line 
metallicity estimates are small (< 0.01 dex), as deter- 
mined by propagating the errors in the line flux mea- 
surement through the abundance calculation. However, 
for example the representative systematic error for the 
PP04 03N2 abundance diagnostic is ~ 0.07 dex, as de- 
termined by IKewlev fc Ellisonl ( 120081 ) via comparison to 
other strong line abundance indicators. 

Fortunately, our detection of the weak [O III] A4363 
auroral line (S/N~ 6; Figure 3) allows us to derive an 
oxygen abundance via the “direct,” T e method. We em- 
ploy ajnethodologv similar to that used by, for exam- 
ple, iLevesauc et all (12010a ). We first derive the elec- 
tron temperature (T e = 1.09 ±0.03 x 10 4 K) and density 
(n e = 83.0 ±7.0 cm -3 ) from the [O III] and [S II] line ra- 
tios using the temden task of the IRAF package nebular 
dShaw fc Dufourlll994f). der i ve the 0 + temperature using 
the calibration of iGarnettl (|1992f l. an d finally e s timat e 
the 0+ and 0 ++ abundances following iShi et all (120061) . 
The direct abundance, log(0/H)±12 =8.24±0.03, is in 
good agreement with the PP04 03N2 value. The stated 
uncertainty reflects the propagation of the uncertainties 
for the line flux measurements. Indeed, the offset be- 
tween thes e two diagnostics should be very small at this 
metallicity IKewlev fc Ellisonl (|2008l ). 

We estimate the star formatio n rate (SF R) of th e host 
galaxy using the Ha relation of iKcnnicuttl (I1998D . Af- 
ter correcting for host galaxy extinction, we measure the 
Ha luminosity from our Gemini spectrum (Table [3]) and 
estimate SFR= 1.1 M 0 yr _1 . 

6.1. Blue Compact Galaxy Host 
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We compare the host galaxy of SN 2010ay to the 
nearby galaxy population of the SDSS spectroscopic sur- 
vey. The physical properties of the host galaxy, SDSS 
J123527. 19+270402. 7, are estimated in the MPA/JHU 
cataloEpl. The total (photometric) galaxy stellar mass 
(M*) is given as 3.61^3 x 10 8 M 0 , the aperture-corrected 
SFR is l-0j^2 Af© yr -1 , and the nuclear (fiber) oxy- 
gen abundance (0/H o ) is log(0/H)+12 =8.58 ^q;q 3 on 
the scale of ITremonti et all (12004T) (T04). The specific 
star formation rate (SSFR) of the host galaxy is then 
ss 2.8 ^q 4 Gyr -1 . For consistency, we consider these val- 
ues of M*, the oxygen abundance, and the SFR for the 
host galaxy of SN 2010ay when comparing to other galax- 
ies in the MPA/JHU catalog. 

The oxygen abundance and SFR of the host galaxy 
of SN 2010ay listed in the MPA/JHU catalog are con- 
siste nt with the values we derive in this paper (see 
also iKellv fe Kirshneil 1201 111 . The MPA/JHU cata- 
log lists metallic i ties o n the T04 scale. Using the 
IKewlev fc Ellisonl (|2008l ) conversion to the PP04 scale, 
the T04 metallicity estimate corresponds to a metallic- 
ity of log(0/H)+12 = 8.38, which is ~ 0.2 dex higher 
than the one we measure (log(0/H)+12 = 8.19). How- 
ever, there is a large (~ 0.2 dex) rms scatter between 
the PP04 03N2 and T0 4 diagnostics at the reg ime of 
log(0/H)+12ppo4 ~ 8.2 (IKewlev fo Ellisonl 120081) . The 
SFR in the MPA/JHU catalog is also in good agreement 
with the value we estimate from the Ha luminosity. Al- 
though out estimate does not include an aperture cor- 
rection, the size of the Gemini slit (1") should encom- 
pass most of the star formation in the galaxy (Petrosian 
r = 1.355"; Tabled]). 

The mass to light ratio of the host galaxy of SN 2010ay 
is low compared to typical star-forming galaxies. To 
compare the host galaxy to the general galaxy popu- 
lation, we select a subset of the MPA/JHU catalog by 
requiring that estimates of Af * , SFR, and O/l I, , be avail- 
able a nd we remove AGN according iKauffmann et all 
(120031) . We consider 167,837 starbursting galaxies fol- 
lowing these constraints. The host galaxy ranks in 
the [4th,38th, 1 ltli] percentile in [M*,SFR,0/H o ] among 
these galaxies. Among the selected galaxies with a stel- 
lar mass as low as the host galax-so, the median and 
standard deviation of the B-banc0 absolute magnitude 
is —15.8 ± 1.3 mag. With Mb =—18.35 ± 0.05!mag, the 
host of SN 2010ay is brighter than other galaxies with a 
similar mass at the 2er level. The discrepancy is due to 
the blue color of SN 2010ay, which indicates a best-fit sin- 
gle stellar population that is very young and therefore has 
a low stellar mass to light ratio. Among the 1,184 galax- 
ies in the MPA/JHU catalog that meet the constraints 
above and have a color similar to the host of SN 2010ay 
(0.47 < u — r < 0.67, from SDSS fiber magnitudes), the 
host galaxy has typical properties, with [M*,SFR,0/H 0 ] 
in the [46th,49th54th] percentile. 

23 http:/ /www.mp a -garching. mpg.de/SDS S (described in 
Kauffmamve^sH 12003; ITremonti et alj 120041 : iBrinchma nn et al.l 
2003 T l^aImi'et*aLl [2(1571 . and updated for SDSS DR7) 

24 This subset is selected such that the host galaxy of SN 2010ay 
has the median mass: 1.39 < M* < 5.62 X 10 8 Mq, 7V s j m = 6, 978. 

25 We obtain B-band magnitudes by converting the /^-corrected 
gri magnitudes given in the MPA/JHU catalog t o BVR , magnitudes 
using the transformation of IBlanton fc Roweisl 1120071 ). 


Based on these properties, we classify the host galaxy 
of SN 2010ay as a luminous Blue Compact Galaxy 
(BCG). BCGs span a large range in luminosity (—21 < 
Mb < —12, where luminous BCGs have Mb < —17), but 
are distinguished by their blue colors (B — V < 0.45), 
high SFR (1 < SFR < 20 M 0 yr -1 ), and low metal- 

licit^ (Z^/50 < Z < Z®/ 2; iKunth fc Ostlinl 120001 : 

iKong fc Chend 120021 ). The host galaxy of SN 2010ay 
has a luminosity (Mb =—18.35 ± 0.05), color (B — V = 
0.11 ± 0.07), SFR (I.O +02 Af© yr -1 ), and metallicity 
(Z ~ 0.3 Z Q ) consistent with all these ranges 0 

6.2. Comparison to SNe Ic-BL and GRB-SNe Host 
Galaxies 

Our measurement of the metallicity from the Gemini 
spectrum indicates that the explosion site of SN 2010ay is 
~ 0.2(0. 5) dex lower in metalli city than the me dian SNe 
Ic (Ic-BL) in the sample of lModiaz et al.l (|2011l) . In that 
sample, the median PP04 03N2 metallicity measured at 
the explosion site of SNe Ic is log(0/H)+12 « 8.7 and for 
Ic-BL is ss 8.4 dex, for 12 and 13 objects, respectively. If 
instead the KD02 metallicity is used, the median of the 
sample is « 8.9 dex for SNe Ic (13 objects) and « 8.7 dex 
for Ic-BL (15 objects), so the abundance of the SN 2010ay 
host galaxy is similarly low compared to the median. 

The metallicity of the environment of SN 2010ay 
is more similar to previously-studied nearby GRB-SN 
progenitors. A metallicity identical to our measure- 
ment was measured a t the explosion site of SN 2010bh 
dLevesaue et al.ll2011l) : log(Q /H)+12 = 8.2. In the sur- 
vev of lLevesaue et all (|2010al) . and adding the measure- 
ment for SN 2010bh, the GRB-SNe host galaxies have an 
average and standard deviation PP04 03N2 metallicity 
of log(0/H)+12 = 8.1 ± 0.1 on the PP04 scale, which is 
consistent with the SN 2010ay env ironment. Amo ng the 
17 LGRB host galaxies surveyed in lSavaglio et al.l (I2009D 
the average metallicity is somewhat lower, 1/6 Z © or 
log(0/H)+12 ~ 7.9, but these are at an average redshift 
of ^ ~ 0.5 that is much higher than SN 2010ay. 

This evidence suggests that the host galaxy of 
SN 2010ay has chemical properties more consistent with 
LGRBs/GRB-SNe than SNe Ic-BL without associated 
GRBs; however, selection effects may mitigate this dis- 
crepancy. SNe found in targeted surveys of bright galax- 
ies have host galaxy properties biased towards higher 
metallici ties, due to the lumin osity-metallicity (L — Z) 
relation (ITremonti et nllfjQO 1 :. LGRBs are found in an 
untargeted manner through their gamma-ray emission 
and therefore are not biased by this relation. 

SN 2010ay joins a growing list of SNe Ic-BL that 
have been discovered in low metallicity host galaxies. 
Given the systematic uncertainty in strong line oxygen 
abundance diagnostics (~ 0.07 dex), we will consider 
host galaxies with metallicity Iog(0/H)ppo4 + 12 < 8.3 
(Z < OAZq) to be in the low-metallicity regime of 
SN 2010ay. Among the 15 SNe Ic-BL (9 discovered by un - 
targ eted searches) i n the surveys of lModiaz et all (|2008l ) 
and iModiaz et al.l (I2011D . 4 were found in low metal- 

26 We note that a large fraction of luminous BCGs show evidence 
for disturbed morph o logies or interaction with close companions 
(I Garland et al.l 120041 ; I L6pez- Sanchez et al.l 120061) , but we do not 
see evidence for a companion at the depth of SDSS images of the 
host galaxy of SN 2010ay. 
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licity environments: SN [2007eb,2007qw,2005kr,2006nx] 
at log(O/H)pp 0 4 + 12 =[8.26,8.19,8.24,8.24], All of 
these SNe w ere d iscovered by untargeted searches. 
I Young et al ., 1 (l201fll ) measure the metallicity of the 
host galaxy of the broad-lined Ic SN 2007bg to be 
log(O/H)pp 0 4 + 12 = 8.18, although this SN has 
lightcurve and spectral proper ties that distinguish 
it from normal SNe Ic-BL ( 113.31) . Furthermore, 
lArcavi et al.1 (I2010T) find that SNe Ic-BL are more com- 
mon in dwarf (M r >= —18) host galaxies, which the 
authors attribute to a preference for lower metallicities. 

The star formation properties of the host galaxy of 
SN 2010ay also resemble the host galaxies of LGRBs. 
If we consider those galaxies in the MPA/JHU catalog 
with masses similar to the host galaxy of SN 2010ay 
(as defined above), then the median SFR and 0/H o of 
these galaxies is 0.13 M 0 yr -1 and log(0/H o )+12 = 
8.36, respectively. The host galaxy of SN 2010ay is in 
the [96th, 77th] percentile for [SFR,0/H o ] among these 
galaxies. This indicates that, while the host galaxy of 
SN 2010ay falls within lcr of the mass-metallicity ( M — Z ) 
relation for star-forming galaxies, its SFR is extreme for 
its mass. The 39LGRB host galaxies in the survey of 
iSavaglio et al.l (l2009f ) are similarly low in mass and have 
high star formation rates, with an average stellar mass 
of M* ~ 10 9 Mq and SSFR- 3.5 Gyr” 1 . 

The host galaxy of SN 2010ay falls below the L — Z 
relation for nearby star-forming galaxies, as illustrated 
by Figure [9] We have transformed the fc-corrected 
nri magnitudes from the MPA/JHU catalog to B-band 
felanton fe Roweisll2007D . At the luminosity of the host 
galaxy of SN 2010ay, the median metallicity and stan- 
dard deviation of the SDSS galaxies on the T04 scale is 
log(0/H)+12 =8.93 ±0.17; the host galaxy of SN 2010ay 
falls in the 3rd percentile. In other words, the host galaxy 
of SN 201 0ay is a 2 a outlier fro m th e L — Z relation . 
Similarly, iLevesaue et al.l (12010a!) and lHan et ah! (12010T ) 
suggest that the host galaxies of LGRBs fall below the 
L—Z relation as defined by normal star-forming galaxies, 
B CGs, and the ho s t gala xies of Type Ic SNe. 

iMannucci et al.l (| 201 If ) have explained the offset of 
LGRB host galaxies from the M — Z relation as a 
preference for LGRBs to occur in host galaxies with 
high SFR, as characterized by the Fundam e ntal M etal- 
licit y Relation (FMR) of Mannucci et all (|2010H (see 
also lLara-Lopez et al.ll2010D . Using the extended FMR 
for low mass galaxies from (IMannucci et al.1 120111) . the 
host galaxy of SN 2010ay should have a metallicity of 
log(0/H)+12 = 8.20 give n its stellar mass an d SFR. The 
FMR is calibrated to the iNagao et all (120061 ) metallicity 
scale, which is similar to that of PP04 at this metallic- 
ity. Given the intrinsic scatter in the extended FMR on 
the order of ~ 0.05 dex, this value is consistent with 
the PP04 value we m easure from our G e mini spectrum: 
log(0/H)±12 = 8.19. iKocevski fc Westl (120111 ) similarly 
explain the offset of LGRB host-galaxies from the M — Z 
relation as a SFR effect, but suggest that the long GRB 
host galaxies have even higher SFR than would be im- 
plied by the FMR. 

SN 2010ay is an example of a SN Ic-BL where the host 
galaxy is consistent with the M — Z relation for star- 
forming galaxies, but deviates from the L — Z relation 
due to its low stellar mass to light ratio (1 16.11) . Its 2cr 


TABLE 5 

Comparison between SN2010ay and GRB 
100316D/SN 2010bh 


Property SN2010bh SN2010ay 


Host galaxy properties 

log(0/H)+12 a 8^2 8.19 

Redshift (z) 0.059 0.06717 

M r -18.5 -18.94 


Explosion properties 


SN type 

Ic-BL 

Ic-BL 

v ph (10 3 s -1 ) b 

24 

22 

Mr 

-18.60 ±0.08 

-20.2 ±0.2 

M Ni ( Mq ) 

0.10 ±0.01 

0 9 +01 
u * y -0.i 

Me j (Mq) 

1.93-2.24 

> 4.7 

Ek,51 

12.0-13.9 

> 10.8 

GRB energy (E iso , erg) 

> 5.9 x 10 49 c 

< 6 x 10 48 


Note. — The observe d properties of SN 2010 bh and its 
host galaxy are given by IChornock et al . ( 120101 ) and light 
curve modelling was performed bv ICano et all J2011a[) . The 
properties of SN 2010ay are derived in this paper. 
a The oxygen abundance derived from the PP04 03N2 
metallicity diagnostic, as discussed in fj6] 
b The photospheric expansion velocity at 30 days after ex- 
plosion, as measured from the Si II A6355A feature in U] 
c The lower limit of the total isotropic energy release esti- 
mated bv IStarling et all i ]20i 'll) . 

discrepancy from the L — Z relation would be hard to 
explain as a SFR rate effect alone because among galax- 
ies in the MPA/JHU catalog without AGN (as defined 
above) and with Mb within 0.1 mag of the host galaxy 
of SN 2010ay, the host galaxy has a SFR in the 26th per- 
centile (< lcr discrepancy). 

7. DISCUSSION 

SN 2010ay has all the hallmark features associated with 
GRB- SNe, and yet we find no evidence of a relativis- 
tic explosion. We are able to place constraints on the 
energy, density, velocity, progenitor mass-loss rate, and 
gamma-ray flux of any GRB that may have been associ- 
ated with it. In particular, we may rule out the associa- 
tion of a GRB that looks similar to any spectroscopically 
confirmed GRB-SN to date, except for XRF 060218. 

The low metallicity of the host environment of 
SN 2010ay is inconsistent with the expectation of 
GR B jet formation i n the “collapsar” model. In 
iMacFadven fc Wooslevl (fl999l ) , a high rate of rotation in 
the core of the progenitor is required to power a relativis- 
tic jet. A low metallicity is prescribed to suppress the line 
driven winds that would deprive the core of angular mo- 
mentum. Previous surveys have found observational evi- 
dence for the preferential occurrence of GRB-SN e in low- 
metallicity host galaxies (ILevesaue et al.lfeOlOa ii). Chal- 
lenging this view is the recent discovery of SN 2009bb, 
a broad-lined, engine-driven Type-Ic supernova found 
in a high-metallicity host environmen t ( Sodc rberg et all 
12010b ILevesaue et al.ll2010cl iPignata et al 11201 ll ). In the 
case of SN 2010ay, we have found the opposite case - a 
broad-lined Type Ic supernova found in a low-metallicity 
host environment, but without any indication (via either 
radio or gamma-ray emission) of a central engine. The 
existence of SNe 2009bb and 2010ay emphasizes that pro- 
genitor metallicity is not the key factor that distinguishes 
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Fig. 9. — A plot of host galaxy metallicity versus absolute B magnitude for SNe Ic-BL (red) and engine-driven explosions (blue). The L-Z 
relation of nearby starforming galaxies is plotted as a solid line, with the 15th and 85th percentile boundaries o f the galaxy distribution 
(dashed lines). Here we have transformed the fc-correc ted gri magnitudes from the MPA/JHU catalog to B-band IIBlanton fc Roweisll2007f ) 
and converted the metallicity values to the PP04 scale IIKewlev Hz EllisorJ2008ll . for the purpose of comparing it to metallicity measurements 
for SN host galaxies in the literature. The dot-dashed line is the divider between GRB-SNe and SNe Ic-BL host env ironments suggested 
bvIMod iaz ct al. (2008). Th e host galaxy prope rties of GRB/SNe other than 20 10ay are from the following reference s: [Starling ct aU I201 II : 
IGano ct al. 2011a] ('2010bh'). IModiaz et al.ll2011l ('Ic-BL'I. ILevesque et al.ll2010dl (2009bb), and lLevesque et al.N2010al (other GRB-SNe). 


GRB-SNe from broad-lined SNe Ic without associated 
relativistic ejecta. 

We compare the photospheric expansion velocity of 
SNe Ic-BL and GRB-SNe to the metallicity of their host 
environments in Figure [TOl This comparison emphasizes 
the diversity of explosion and host galaxy properties ob- 
served in both GRB-SNe and SNe Ic-BL. The GRB-SNe 
with the largest velocity gradients (2003dh, a = —0.9; 
2009bb, a = —0.9) occur at metallicities different by 
an order of magnitude (2003dh, log(0/H)+12 = 8.0; 
2009bb, log(0/H)+12 = 9.0). Furthermore, the ve- 
locity gradient of GRB-SNe seems to be uncorrelated 
with the velocity at early times; 2010bh ( a = —0.2) 
and 2006aj (a = —0.3) have similar velocity gradients, 
but 2010bh was ~ 10, 000 km s -1 faster at 10 days 
after explosion. Among SNe Ic-BL from low metallic- 
ity environments (log(0/H)+12 < 8.5), there is a large 
range in both the characteristic velocity and velocity 
gradient ( v ^ = [10, 7, 22] and a = [—0.5, —0.2, —0.4] 
for SNe [2003jd,2007bg,2010ay]). For the three SNe Ic- 
BL at higher metallicities, the characteristic velocities 
tend to be slower (but more consistent) and the veloc- 
ity gradients are stronger are more tightly distributed 
( v ph = [6,4,10] and a = [—0.8, —1.1, —0.5] for SNe 
[1997ef,2002ap,2007ru]). However, a larger sample is 


needed to exclude the possibility of SNe Ic-BL from 
super-solar metallicity environments that have high char- 
acteristic velocities or shallow velocity gradients. 

The fact that a GRB was not detected in association 
with SN 2010ay, despite its similarities to the GRB-SNe, 
could indicate that the relativistic jet produced by this 
explosion was “suf focated” before it em erged from the 
progenitor star dMacFadven et al] 120011 ) . In this sce- 
nario, the duration of the accretion event onto the newly- 
formed central engine is short-lived and the jet post- 
breakout outflow is not ultra-relativistic. In the process 
of being suffocated, the jet transfers momentum to the 
ejecta such that the spectrum is broad lined and the pho- 
tospheric velocity is very high, even at late times, as we 
observe (33]) - However, the low host environment metal- 
licity we measure for SN 2010ay, which is similar to GRB- 
SNe, does not suggest high angular momentum loss that 
would help to weaken the jet. Another alternative is that 
binary interaction plays a key role in the commonality of 
relativistic outflows among SNe Ic-BL. 

Looking forward, additional SNe Ic-BL in sub-solar 
metallicity host environments will be found through 
current and future generations of untargeted transient 
searches. Untargeted searches are required to find su- 
pernovae from low-metallicity host environments, be- 
cause targeted searches only probe the highest metallic- 
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Fig. 10. — SN photspheric velocity, as traced by the Si II A6355A feature, versus host galaxy oxygen abundance for SN 2010ay and 
other Ic-BL (red) and SNe-GRB (blue) from the literature (as in Figure |5|) . The range of velocities hatched for each object comes from 
the velocity at 10 days and at 30 days after explosion, according to the powe r law fits pre sented in F i gure | 5l The oxyge n abundance 
measurements using the PP04 03N2 di agnos tic are from I Levesque et al.l (I2010al) f GRB-SNe j JSahu et al.l }2009) (SN2007ru) , I Young et al.l 
(2010) (SN 2007bg), and Modi az ct al. ( 2011) (other SNe Ic-BL). The range of oxygen abundanc e hatched reflects the er ror bars quoted in 
the literature (when stated) plus the ~ 0.07 dex systematic error of the PP04 03N2 diagnostic (|Kewlev fc Ellisonl 1 2008l b 


ity galaxies due to the luminosity-metallic ity relationship 
rtModiaz et alJl20Ffl : iLeloudas et, al.lfTTTTIl . 

8. CONCLUSIONS 

The optical photometric, spectroscopic, and radio ob- 
servations of SN 2010ay presented here provide an exam- 
ple of a Type Ic-BL SN with explosion and host prop- 
erties similar to the known GRB-SN SN 2010bh. This 
object demonstrates that SNe in low-metallicity envi- 
ronments with high-velocity ejecta are not necessarily 
accompanied by the long-lived relativistic jets necessary 
for a GRB. The existence of SN 2010ay and SN 2009bb (a 
central-engine driven event from a high-metallicity host 
environment) indicates that progenitor metallicity may 
not be the key factor that distinguishes GRB-SNe from 
normal broad-lined Type Ic supernovae. 

We conclude that: 


1. Pre-discovery imaging of the SN2010ay from the 
Pan-STARRSl 37r survey allows us to tightly 
constrain the early-time light-curve of SN 2010ay 
(see Figured]) and explosion date (2010 February 
21.2±1.5), allowing us to search for gamma-ray 
emission that may have been associated with the 
explosion. By fittin g the t emplate SN Ibc light- 
curve of lDrout et al.l (j2010n . we derive an R-band 
peak absolute magnitude of — 20.2±0.2 mag - mak- 
ing SN 2010ay among the brightest SNe Ibc ever 
observed. This peak magnitude suggests that a 
large mass of nickel, M m ~ 0.9t 0 M 0 , has been 
synthesized. The ratio of Mm to M e j is > 2 x larger 
than in known GRB-SNe. 

2. Spectroscopy (see Figure [3]) at the explosion site 
in the host galaxy of SN 2010ay indicates that the 
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host environment of the progenitor star had a sig- 
nificantly sub-solar metallicity ( Z ~ 0.3 Zq), sim- 
ilar to the host environments of known GRB-SNe 
progenitors. 

3. The Ic-BL SN 2010ay strongly resembles the GRB- 
SN 2010bh, particularly in light of its unusually 
high photosphcric expansion velocities at late times 
( v p h ~19.2xl0 3 km s _1 at 14 days after peak) and 
low metallicity host environment. The comparison 
between these two SNe is summarized in Table [5] 

4. Non-detections in late-time radio observations of 
the SN with the EVLA rule out the association 
of a GRB of the nature of the spectroscopically- 
confirmed GRB-SNe, except for the radio afterglow 
associated with XRF 060218. Our radio obser- 
vations imply limits on the velocity, energy, and 
density of any associated relativistic jet and the 
mass loss rate of the progenitor (see Figures 01 
1313 and ©■ Additionally, no coincident gamma- 
ray emission was detected by satellites: the non- 
detection by the interplanetary network indicates 
E-y < 6 x 10 48 erg, while the non-detection by 
the Swift BAT indicates that the peak energy of 
the burst was < 1 x 10 47 erg s _1 if the burst oc- 
curred during the ~ 20% of the explosion window 
when it was in the field of view of the instrument. 
This rules out associated prompt emission similar 
to that of GRBs 031203, 030329, or 100316D, but 
not GRBs 980425 or 060218. 

The pre-discovery imaging of SN 2010ay demonstrates 
the capability of the untargeted PS1 survey for identi- 
fying and monitoring exotic transients, not only in its 
high-cadence Medium-Deep Fields, but also in the all-sky 
37 t survey. Additional detections and multi-wavelength 
follow-up observations of SNe Ic-BL will help to illumi- 
nate the role that ejecta velocity and progenitor metal- 
licity play in the GRB-SNe connection. 
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